Cross-Ancestry Genome-Wide Association Study Defines the Extended CYP2D6 Locus as the Principal Genetic Determinant of Endoxifen Plasma Concentrations.

The therapeutic efficacy of tamoxifen is predominantly mediated by its active metabolites 4-hydroxy-tamoxifen and endoxifen, whose formation is catalyzed by the polymorphic cytochrome P450 2D6 (CYP2D6). Yet, known CYP2D6 polymorphisms only partially determine metabolite concentrations in vivo. We performed the first cross-ancestry genome-wide association study with well-characterized patients of European, Middle-Eastern, and Asian descent (n = 497) to identify genetic factors impacting active and parent metabolite formation. Genome-wide significant variants were functionally evaluated in an independent liver cohort (n = 149) and in silico. Metabolite prediction models were validated in two independent European breast cancer cohorts (n = 287, n = 189). Within a single 1-megabase (Mb) region of chromosome 22q13 encompassing the CYP2D6 gene, 589 variants were significantly associated with tamoxifen metabolite concentrations, particularly endoxifen and metabolic ratio (MR) endoxifen/N-desmethyltamoxifen (minimal P = 5.4E-35 and 2.5E-65, respectively). Previously suggested other loci were not confirmed. Functional analyses revealed 66% of associated, mostly intergenic variants to be significantly correlated with hepatic CYP2D6 activity or expression (ρ = 0.35 to -0.52), and six hotspot regions in the extended 22q13 locus impacting gene regulatory function. Machine learning models based on hotspot variants (n = 12) plus CYP2D6 activity score (AS) increased the explained variability (~ 9%) compared with AS alone, explaining up to 49% (median R2 ) and 72% of the variability in endoxifen and MR endoxifen/N-desmethyltamoxifen, respectively. Our findings suggest that the extended CYP2D6 locus at 22q13 is the principal genetic determinant of endoxifen plasma concentration. Long-distance haplotypes connecting CYP2D6 with adjacent regulatory sites and nongenetic factors may account for the unexplained portion of variability.

Pharmacogenetics of UGT1A4, UGT2B7 and UGT2B15 and Their Influence on Tamoxifen Disposition in Asian Breast Cancer Patients.

Tamoxifen (TAM) is an established endocrine treatment for all stages of oestrogen receptor (ER)-positive breast cancer. Its complex metabolism leads to the formation of multiple active and inactive metabolites. One of the main detoxification and elimination pathways of tamoxifen and its active metabolites, 4-hydroxytamoxifen (4-OHT) and endoxifen, is via glucuronidation catalysed by uridine 5'-diphospho-glucuronosyltransferases (UGTs). However, few studies have comprehensively examined the impact of variations in the genes encoding the major hepatic UGTs on the disposition of tamoxifen and its metabolites. In the present study, we systematically sequenced exons, exon/intron boundaries, and flanking regions of UGT1A4, UGT2B7 and UGT2B15 in 240 healthy subjects of different Asian ethnicities (Chinese, Malays and Indians) to identify haplotype tagging single nucleotide polymorphisms. Subsequently, 202 Asian breast cancer patients receiving tamoxifen were genotyped for 50 selected variants in the three UGT genes to comprehensively investigate their associations with steady-state plasma levels of tamoxifen, its active metabolites and their conjugated counterparts. The UGT1A4 haplotype (containing variant 142T>G, L48 V defining the *3 allele) was strongly associated with higher plasma levels of TAM-N-glucuronide, with a twofold higher metabolic ratio of TAM-N-glucuronide/TAM observed in carriers of this haplotype upon covariate adjustment (P < 0.0001). Variants in UGT2B7 were not associated with altered O-glucuronidation of both 4-OHT and endoxifen, while UGT2B15 haplotypes had a modest effect on (E)-endoxifen plasma levels after adjustment for CYP2D6 genotypes. Our findings highlight the influence of UGT1A4 haplotypes on tamoxifen disposition in Asian breast cancer patients, while genetic variants in UGT2B7 and UGT2B15 appear to be of minor importance.

Association of CYP2C19*2 and associated haplotypes with lower norendoxifen concentrations in tamoxifen-treated Asian breast cancer patients.

AIM: The aim was to examine the influence of CYP2C19 variants and associated haplotypes on the disposition of tamoxifen and its metabolites, particularly norendoxifen (NorEND), in Asian patients with breast cancer. METHODS: Sixty-six CYP2C19 polymorphisms were identified in healthy Asians (n = 240), of which 14 were found to be tightly linked with CYP2C19*2, CYP2C19*3 and CYP2C19*17. These 17 SNPs were further genotyped in Asian breast cancer patients receiving tamoxifen (n = 201). Steady-state concentrations of tamoxifen and its metabolites were quantified using liquid chromatography­mass spectrometry. Non-parametric tests and regression methods were implemented to evaluate genotypic­phenotypic associations and haplotypic effects of the SNPs. RESULTS: CYP2C19 functional polymorphisms and their linked SNPs were not significantly associated with plasma concentrations of tamoxifen and its main metabolites N-desmethyltamoxifen, (Z)-4-hydroxytamoxifen and (Z)-Endoxifen. However, CYP2C19*2 and its seven linked SNPs were significantly associated with lower NorEND concentrations, MRNorEND/NDDM and MRNorEND/(Z)-END. Specifically, patients carrying the CYP2C19*2 variant allele A had significantly lower NorEND concentrations [median (range), GG vs. GA vs. AA: 1.51 (0.38­3.28) vs. 1.28 (0.30­3.36) vs. 1.15 ng ml−1 (0.26­2.45, P = 0.010)] as well as significantly lower MRNorEND/(Z)-END [GG vs. GA vs. AA: 9.40 (3.27­28.35) vs. 8.15 (2.67­18.9) vs. 6.06 (4.47­14.6), P < 0.0001] and MRNorEND/NDDM [GG vs. GA vs. AA: 2.75 (0.62­6.26) vs. 2.43 (0.96­4.18) vs. 1.75 (1.10­2.49), P < 0.00001]. CYP2C19 H2 haplotype, which included CYP2C19*2, was also significantly associated with lower NorEND concentrations (P = 0.0020), MRNorEND/NDDM (P < 0.0001) and MRNorEND/(Z)-END (P < 0.0001), indicating significantly lower formation rates of NorEND. CONCLUSION: These data highlight the potential relevance of CYP2C19 pharmacogenetics in influencing NorEND concentrations in tamoxifen-treated patients, which may influence treatment outcomes.

Genetic variations of NR1I3 and NR2B1 in Asian populations.

Several nuclear receptors are being increasingly recognized for their role as master xenosensors. Among them, CAR-RXRα heterodimer, as encoded by NR1I3 and NR2B1, responds to the presence of drug compounds and regulates the transcription of a wide array of genes involved in their disposition. To investigate the frequency distribution and linkage disequilibrium patterns of NR1I3 and NR2B1 genetic variations, these genes were screened in 168 healthy local Asian subjects, namely Chinese, Malays, and Indians (n=56 subjects each). A total of 38 and 88 SNPs were identified in NR1I3 and NR2B1, respectively. Among them, there were 13 and 43 novel SNPs present at low allelic frequencies (<10%) in NR1I3 and NR2B1, respectively. Notably, the genetic variations in the NR1I3 and NR2B1 genes were mainly confined to the introns whilst the exons were highly conserved across the ethnic populations. Indians harboured distinct frequency distributions from Chinese and Malays in both genes. Based on the linkage disequilibrium patterns of both genes, a number of tag-SNPs were selected for each population (n=8-13 for NR1I3; n=12-18 for NR2B1). In-silico prediction analyses revealed a number of possible functional SNPs. Our data would be valuable for future pharmacogenetic studies on the drug substrates of CAR-RXRα target genes.

Genetic variations of NR2A1 in Asian populations: implications in pharmacogenetics studies.

HNF4α (encoded by gene NR2A1) is a dominant transcriptional regulator of various drug disposition genes. It forms a circuitry of molecular cross-talk with other nuclear receptors such as PXR and CAR to synergistically initiate transcription. This study reports on the frequency, linkage disequilibrium pattern and tag-SNP selection of NR2A1 polymorphisms in three local Asian populations, namely Chinese, Malays and Indians (n = 56 subjects each). A total of 69 polymorphisms were identified in the genomic region of NR2A1, of which thirty-three were novel polymorphisms with low allelic frequencies (<0.02). The exonic region of NR2A1 was highly conserved with only 4 novel and 1 reported SNPs identified at low allelic frequencies of less than 0.02. Based on the criteria of MAF ≥ 0.05 and R(2) ≥ 0.80, there were 19, 20 and 22 tag-SNPs selected to represent the genetic polymorphisms of NR2A1 in Chinese, Malays and Indians, respectively. In-silico predictions suggested that some of these polymorphic variants may exert functional effects through affecting the binding sites of transcription and splicing factors. Our study provides valuable information on the genetic variability of NR2A1 which would be useful for pharmacogenetics studies in the local Asian populations.

Pharmacogenetics of CYP1A2, novel polymorphisms and haplotypes in three distinct Asian populations.

CYP1A2 play an important role in the metabolism of many carcinogens and clinically important drugs. CYP1A2 activity has been found to be influenced by the presence of polymorphic variants which were reported to display wide interethnic variation. This study investigates the frequency distribution and linkage disequilibrium patterns of CYP1A2 genetic polymorphisms, and characterize their haplotype structures in three healthy Asian populations in Singapore (Chinese, Malay, and Indian). The entire CYP1A2 gene was screened in 126 healthy subjects from all three ethnic groups (N=42 each). A total of 25 polymorphisms was identified, of which nine were novel. The polymorphisms, -2467delT and -163C>A were detected at high frequencies in all Asian ethnic groups. Significant interethnic differences were observed in the genotypic frequency distribution of IVS2-99G>A (P<0.01) and 1548C>T (P=0.05) across the three ethnic groups while -163C>A (P=0.02) was found to differ between Chinese and Malays. Haplotype analyses revealed four to six major haplotypes in each ethnic population which accounted for more than 60% of the cumulative haplotype frequencies. Future studies should be done to investigate the functional roles of these haplotypes.